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Intrinsic Parameters of GRB990123 from Its Prompt 
Optical Flash and Afterglow 
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ABSTRACT 

We have constrained the intrinsic parameters, such as the magnetic energy density 
fraction (es), the electron energy density fraction (e e ), the initial Lorentz factor (Fo) 
and the Lorentz factor of the reverse external shock (T rs ), of GRB990123, in terms 
of the afterglow information (forward shock model) and the optical flash information 
(reverse shock model). Our result shows: 1) the inferred values of e e and €b are consis- 
tent with the suggestion that they may be universal parameters, comparing to those 
inferred for GRB970508; 2) the reverse external shock may have become relativistic 
before it passed through the ejecta shell. Other instrinsic parameters of GRB990123, 
such as energy contained in the forward shock E and the ambient density n are also 
determined and discussed in this paper. 
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1 INTRODUCTION 



The current standard model for gamma-ray bursts (GRBs) 
and their afterglows is the fireball-plus-shock model (see Pi- 
ran 1999 for a review). It involves that a large amount of en- 
ergy, Eo ~ 10 51-54 ergs, is released within a few seconds in a 
j^ ■ small volume with negligible baryonic load, Mc 2 <C Eo- This 
J— i leads to a fireball that expands ultra-relativistically with a 
Lorentz factor To — Eo/Mc 2 > 100 required to avoid the at- 
tenuation of hard 7-rays due to pair production (e.g. Woods 
& Loeb 1995; Fenimore, Epstein & Ho 1993). A substantial 
fraction of the kinetic energy of the baryons is transferred 
to a non-thermal population of relativistic electrons through 
Fermi acceleration in the shock (Meszaros & Rees 1993) . The 
accelerated electrons cool via synchrotron emission and in- 
verse Compton scattering in the post-shock magnetic fields 
and produce the radiation observed in GRBs and their af- 
terglows (e.g. Katz 1994; Sari et al. 1996; Vietri 1997; Wax- 
man 1997a; Wijers et al. 1997). The shock could be either 
internal due to collisions between fireball shells caused by 
outflow variability (Paczyriski & Xu 1994; Rees & Meszaros 
1994), or external due to the interaction of the fireball with 
the surrounding interstellar medium (ISM; Meszaros & Rees 
1993). The radiation from internal shocks can explain the 
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spectra (Pilla & Loeb 1998) and the fast irregular variability 
of GRBs (Sari & Piran 1997), while the synchrotron emis- 
sion from the external shocks provides a successful model for 
the broken power law spectra and the power law decay of 
afterglow light curves (e.g. Waxman 1997a, b; Wijers, Rees 
& Meszaros 1997; Vietri 1997; Dai & Lu 1998a,b,c; Dai et al. 
1999; Wang et al. 1999a,b; Huang et al. 1998a,b; 1999a,b). 
The properties of the synchrotron emission from GRB 
shocks are determined by the magnetic field strength, B, and 
the electron energy distribution behind the shock. Both of 
them are difficult to estimate from first principles, and so the 
following dimensionless parameters are often used to incor- 
porate modeling uncertainties (Sari et al. 1996), es = —B-, 

e e = — . Here Ub and U e are the magnetic and electron en- 

e th 

ergy densities and eth = nm v c 2 (^ v — 1) is the total thermal 
energy density behind the shocks, where m p is the proton 
mass, n is the proton number density, and y p is the mean 
thermal Lorentz factor of the protons. In spite of these un- 
certainties, an important assumption that es and e e do not 
change with time, has been made in the standard external 
shock model. Through the computation, Wijers & Galama 
(1998; hereafter WG99) even suggested that they may be 
universal parameters, i.e. the same for different bursts. 

The BeppoSAX satellite ushered in 1999 with the dis- 
covery of GRB990123 (Heise et al. 1999), the brightest GRB 
seen by BeppoSAX to date. This is a very strong burst. 
An assumption of isotropic emission and the detection of 
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the source's redshift z = 1.6004 , lead to a huge energy 
release about 1.6 x 10 54 ergs (Briggs et al. 1999; Kulkarni 
ct al. 1999a) in 7-rays alone. GRB990123 would have been 
amongest the most exciting GRBs even just on the basis of 
these facts. Furthermore, ROTSE discovered a prompt opti- 
cal flash of 9-th magnitude (Akerlof et al. 1999). It is the first 
time that a prompt emission in another wavelength apart 
from 7-rays has been detected from GRB. Such a strong 
optical flash was predicted to arise from a reverse external 
shock propagating into the relativistic ejecta (Meszaros & 
Rees 1997; Sari & Piran 1999a,b, hereafter SP99a,b). This 
is the so called "early afterglow" . The five last exposures of 
ROTSE show a power law decay with a slope of ~2.0, which 
can also be explained by the reverse shock model (SP99b; 
Meszaros & Rees 1999). The usual afterglows in X-ray, op- 
tical, IR and radio bands were also detected after the burst. 
They have two distinguishing features: 1) the radio emission 
is unique both due to its very early appearance and its rapid 
decline; 2) the temporal decaying behaviour of the r-band 
light curve after two days steepens from about t~ 1A to £~ 1,8 
(Kulkarni et al. 1999a; Fruchter et al. 1999; Casrto-Tirado 
ct al. 1999), and this steepening might be due to a jet which 
has transited from a spherical-like phase to sideways expan- 
sion phase (Rhoads 1999; Sari et al. 1999; Wei & Lu 1999) 
or a dense medium which has slowed down the relativistic 
expansion of shock quickly to a non-relativistic one (Dai & 
Lu 1999a,b). 

Galama et al. (1999) assumed that the radio emission 
is produced by the forward shock as usual and then recon- 
structed the radio-to- X-ray afterglow spectrum on January 
24.65 UT. However, later work shows that the simplest inter- 
pretation of this "radio flare" is that it arises in the reverse 
shock and that such radio emission is an inevitable conse- 
quence of the prompt bright optical flash seen by ROTSE 
(Kulkarni et al. 1999b; Sari & Piran 1999b). Kulkarni et 
al. (1999b) also constrained two key parameters of the for- 
ward shock, the peak flux F Vm and the peak frequency i> m , 
to within a factor of two. For previous bursts, we have no 
other information apart from the afterglow to infer the in- 
trinsic parameters of the external shocks. But now the opti- 
cal flash of GRB990123 has been fortunately detected, which 
enables us to determine another two key parameters, the 
initial Lorentz factor To and the Lorentz factor of the re- 
verse shock r r3 . WG99 computed the intrinsic parameters 
of GRB970508 and GRB971214 in terms of their afterglow 
data, and found that e e is nearly the same for these two 
bursts, suggesting it may be a universal parameter. Gra- 
not, Piran & Sari (1998a,b; hereafter GPS99a,b) modified 
the set of equations derived by WG99, and inferred the elec- 
tron energy density fraction and the magnetic energy density 
fraction of GRB970508 to be: e e = 0.57; e B = 8.2 x 10~ 3 . 
Here we apply the set of equations of GPS99a to GRB990123 
and try to determine some intrinsic parameters based on the 
information of the two aspects of GRB990123 — the optical 
flash and the afterglow. 

The initial Lorentz factor To is also an important phys- 
ical parameter of GRBs. It is a crucial ingredient for con- 
straining models of the source itself, since it specifies how 
"clean" the fireball is as the baryonic load is M ~ Eq/Tqc 2 . 



Unfortunately, the spectrum of GRBs can provide only a 
lower limit to this Lorentz factor (Fo > 100). Moreover, the 
current afterglow observations, which detect radiation from 
several hours after the burst, do not provide a verification 
of the initial extreme relativistic motion. A possible method 
to estimate To of GRBs has been suggested by Sari & Pi- 
ran (1999a), based on identifying the "early afterglow" peak 
time. In this paper, the initial Lorentz factor has been in- 
ferred more precisely from the full set of equations describing 
the reverse shock region. 

In section 2 and 3, we compute the intrinsic parameters 
of GRB990123 from its afterglow and optical flash informa- 
tion. In the final section, we give our conclusions. 



2 PARAMETERS FROM THE AFTERGLOW 

Intrinsic parameters like the magnetic energy density frac- 
tion tB, electron energy density fraction e e , energy in the 
forward external shock E = E52 x 10 52 ergs and ambient 
density n can be determined from the afterglow spectrum 
(GPS99a,b; WG99), i. e. if we know all three break frequen- 
cies (not necessary at the same time) and the peak flux of 
the afterglow, we can infer all these parameters. 

^From the observations of the afterglow of GRB990123, 
Kulkarni et al. (1999b) have estimated two key parameters 
of the forward shock: u m ~ 1.1 x 10 13 Hz and F Vm ~ 170 /iJy 
at the time t* = 1.25 days after the burst. The cooling fre- 
quency v c cannot be seen from the radio-to-X-ray spectrum 
obtained by Galama et al (1999). This indicates that v c is 
at or above the X-ray frequencies. We need to determine it 
more precisely. The X-ray afterglow, observed 6 hours after 
the burst, decayed with a x = 1.44±0.07 ( Heise et al. 1999), 
while the optical afterglow with a r = 1.10 ± 0.03 (Kulkarni 
et al., 1999). An X-ray afterglow decay slope steeper by | 
than an optical decay, which seems to be the case in this 
burst, is predicted by Sari, Piran and Narayan (1998), if 
the cooling frequency is between the X-rays and the opti- 
cal. So at the time 6 hours after the burst , 4 x 10 14 < v c < 
fx ~ (4.4 — 44) x 10 17 Hz. Extrapolating it to the time i* as 
v c oc £~3, we get 1.8 x 10 14 Hz< v c (U) < (2-20) x 10 17 Hz. 
Another speculative constraint on v c is obtained from the 
GRB spectrum itself by Sari & Piran (1999b), who con- 
strained u c > 2x 10 1 Hz at the time t ~ 50 sec. Extrapolat- 
ing it to t*, we get v c (t*) > 0.4x 10 18 . Now v c is almost deter- 
mined, and we take the approximate value v c = 0.5 x 10 18 Hz, 
which is in agreement with the estimate of Galama et al 
(1999). In addition, Kulkarni et al. (1999a) have inferred the 
electron index p (defined as N(-y e ) oc 7 e ~ p ) to be p = 2.44. 
Now, apart from the self-absorption frequency v a , we have 
all other three quantities of the afterglow spectrum required 
to calculate the intrinsic parameters: 



v m ~ l.lxl0 13 Hz, ^ c ~0.5xl0 18 Hz, F„„ 



170 ^Jy, p~ 2.44(1) 



Following GPS99b, we adopt the formulas of v m and 
Fv m f rom GPS99a and v c from Sari, Piran & Narayan 
(1998): 

~ on^iniS/i 1 \l/2/P~2 2 1 , 2 2,-, 1/2, -3/2tt , n \ 

z/ m ~ 2.9x10 (1+2) ( -) e B e e E 5 2 t d Hz;(2) 

p- 1 
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F Vm ~ 1.7 x 10 4 (1 + z)e B 1 / 2 E 52 ni 1/2 (-^ ^ 2 

10 28 c 



)"Vy; (3) 



m c ~ 2.7 x 10 12 e s - 3/2 £ 5 2- 1/2 n- 1 t d - 1/2 (l + z)- 1/2 Kz, (4) 

where z is the redshift of the burst and cLl — 2c(l + 2: — 
a/1 + z)/Ho is the luminosity distance. 

Now we have three equations (2), (3) and (4) with four 
unknowns: e e , es, -E52 and n. To solve these equations, we 
here assume that the electron density fraction e e is the same 
for different bursts, just an argument of WG99, though in 
which a different set of equations are used. Since we here use 
the formulas of GPS99a, we adopt the value of e e from that 
of GRB970508 inferred according to the above formulas, i.e. 
e e ~ 0.57(GPS99b). 

By combining Eqs.(l) and (2)-(4), we get the values of 
four intrinsic parameters of the forward shock region: 



0.57 t B ~ 3.1 x 10" 



0.01 



(5) 



We astonishingly find that the value e B inferred is very 
close to that of GRB970508 (e s ^ 8.2 x 10" 3 ; GPS99b), 
considering the uncertainties in the v m and F Vm of a factor 
of two. This result supports our above adoption of the value 
of e e and implies that e e and e B may be universal parameters 
(i.e. constants for every GRBs), favouring the argument of 
WG99. 

The ambient density n inferred for GRB990123 is n ~ 
0.01. This density is on the low side of normal for a disc 
of galaxy but definitely higher than expected for a halo, 
lending a further support to the notion that bursts occur 
in gas-rich environment. The inferred isotropic energy left 
in the adiabatic forward shock is E ~ 5 x 10 52 ergs, about 
thirty times less than the isotropic energy in 7-rays. This 
case is very similar to GRB971214 (WG99). We think, for 
GRB990123, there are two processes causing E52 < £52,7- 
One is that there is a radiative evolution phase before the 
adiabatic phase, causing it to emit more of the initial ex- 
plosion energy and leaving less for the adiabatic phase. Ac- 
cording to Sari, Piran & Narayan (1998), we can estimate 
the reduced energy of the forward shock in the self-similar 
deceleration stage, i.e. 



E 



/,52 



0.02e B - 3/5 e e - [i/5 E l , 52 4/5 (r A /W0)- 4/T, n- 2/r >, (6) 



where £^52 denotes the initial isotropic energy of the for- 
ward shock during the self-similar stage and -E/,52 denotes 
the final energy after the radiative phase. The value of £4,52 
is difficult to be determined precisely and it may be less than 
£52,7 for GRB990123, according to Freedman & Waxman 
(2000). If we use £;, 52 ~ £52,7 , e e ~ 0.57, e B ~ 8.2 x 10~ 3 
(instead of e B ~ 3.1 x 10~ 3 , since the latter is less reliable 
than the former, considering the uncertainties in the v m and 
F Um of a factor of two for GRB990123.), Y A ~ 300 (see the 
next section) and n ~ 0.01, we then get -E/,52 ~ 56. The 
other process is the sideways expansion of the fireball jet 
(Kulkarni et al.l999a,b), which can also reduce the energy 
per solid angle, hence the isotropic energy £52. Since the 



opening angle of the jet is 



~r Cgtp 



r/ct 



% + 7~ 



(Sari, Piran & Halpern 1999; Rhoads 1999), at the time 
t = 1.25 days (very near the break time of the jet evolution 
t b ~ 2.1 days), 7 ~ ±- (note that 7 oc i~ 3/8 ), then 9 ~ 2$ . 



Therefore, the real isotropic energy £52 left in the late adi- 
abatic forward shock should be 56(-g- 2 -)~ 2 ~ 14 , in rough 
agreement with the above value inferred from the afterglow 
spectrum. An additional possible loss of energy may be the 
reverse shock itself if it is radiative. 



3 PARAMETERS FROM THE OPTICAL 
FLASH 

An optical flash is considered to be produced by the re- 
verse external shock, which heats up the shell's matter and 
accelerates its electrons (SP99b; Meszaros & Rees 1999). 
BATSE's observations triggered ROTSE via BACODINE 
system (Akerlof et al. 1999). An 11.82 magnitude optical 
flash was detected on the first 5 seconds exposure, 22.18 
seconds after the onset of the burst. Then the optical emis- 
sion peaked in the following 5 seconds exposure, 25 seconds 
later, which revealed an 8.95 magnitude signal (~Uy). The 
optical signal decayed to 10.08 magnitude 25 seconds later 
and continued to decay down to 14.53 magnitude in the 
subsequent three 75 seconds exposures that took place up 
to 10 minutes after the burst. The five last exposures depict 
a power law decay with a slope ~ 2.0 (Akerlof et al. 1999; 
SP99b). Sari & Piran (1999b) and Meszdros & Rees (1999) 
assumed that the ejecta shell follows the Blandford-McKee 
(1976) self-similar solution after the reverse shock has passed 
through it and explained the decay of i~ 2,0 . 

So we assume that at the optical emission peak time 
(t = 50 sec) the reverse shock had just passed through the 
ejecta shell. At this time, the Lorentz factor of the reverse 
shock P rs is approximately given by 



1° 

T A ' 



(7) 



where Po is the initial Lorentz factor of the ejecta and Ta 
is the Lorentz factor of the ejecta at the optical flash peak 
time. Then the random minimum Lorentz factor jmin of the 
electrons in the reverse shock region is: 



^fmin 



m p p- 



1°. 

Ta : 



, . (8) 

m e p- 

The formulas of u m at the reverse shock were given by 
Sari & Piran (1999b). We here add the correction for red- 
shift. 



1.2xl0-(^) 2 (l ^ 



^1/2, £0,2 1/2 

' l 300 j ' 



l+z) _1 Hz< 5xl0 14 Hz.(9) 



The observed flux at v m can be obtained by assuming 
that all the electrons in the reverse shock region contribute 
the same average power per unit frequency P' Vm at v m , which 
is given by P' Vm = — f - , where B' = r 'AC^?,2imm p e B . 
Adding one factor of Ya to transform to the observer frame 
and accounting for the redshift, we have: 

N e V A Pl m {l + z) 



F Vm = 



(10) 



where N e is the total number of radiating electrons in the 
ejecta shell, and dh = 2c(l + z — \f\ + z)/Ho is the lumi- 
nosity distance. Please note that N e here is different from 
the N e adopted in the forward shock region, which is the 
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total number of swept-up electrons by the forward exter- 
nal shock. We consider N e here to be the total number of 
electrons contained in the baryonic load: 



N e = 



Al 



E': 



T m p c 2 



1.08xl0 54 ( 



x-l 



1.6 x 10 5 



3iT>(inSn> .(") 



1000 



where E 1 is the total energy in 7-rays. Substituting the ex- 
pression of P' Vm and 2 = 1.6 into Eq.(10), we get 



F„ 



0.74( 



E- f .. To ._i / .r j 4. 2 1/2/ £S \1/2 T /-iox 

I^To^W W n ( IF^ Jy ' (12) 



Please note that this formula always holds whether the 
ejecta is jet-like or spherical, because the beaming factor 
in Eq.(10) and (11) will cancel out each other in the jet- 
like case. According to the jump condition of the shock, the 
Lorentz factor of the shocked shell should be approximately 
equal to that of the shocked ISM (Piran 1999). The Lorentz 
factor of the forward shocked ISM can be obtained from 
the standard afterglow model (e. g. Sari, Piran & Narayan 
1998): 



T^ /t\ ~ c/-^ 52 \l/8r td -,-3/8 

r AJS (t)^6(—)<(—) 

For £52 ~ 5 and n ~ 0.01, we get 
r A ,r S (50s) =r A ,/ s (50s) =300 



(13) 



(14) 



Taking the above inferred value es ~ 3.1 x 10 3 , from the 
Eqs.(9), (12), and (14) with the conditions: v m < 5 x 10 14 Hz 
and F Vm ~ Uy, we finally get: 



T A = 300, To = 1200, e e < 0.6. 



(15) 



Please note that here the value e e < 0.6 inferred from 
the optical flash data is consistent with that inferred in- 
dependently from the afterglow information. On the other 
hand, if we substitute the value e e — 0.57 into Eq.(9), we find 
that the peak frequency of the reverse shock v m is almost 
located at the optical band. In addition, our inferred initial 
Lorentz factor To is six times larger than that obtained by 
Sari & Piran (1999b), who have used the ambient density n 
of GRB970508. Consequently, at the time the reverse shock 
has just passed through the ejecta shell, its Lorentz factor 
was r rs = -pr 2 - ~ 4. This indicates that the reverse shock had 

1 A 

become relativistic before it crossed the entire shell. This 
result is also different from that obtained by Sari & Piran 
(1999b), who found the Lorentz factor of the reverse shock of 
GRB990123 was only near one. However we argue that our 
result is reasonable according to the criterion presented by 
Sari & Piran (1995) (also see Kobayashi et al. 1998). They 
defined a dimensionless parameter £ constructed from I, A 
and To: 



£ = (VA) 1/2 r - 4/3 , 

where I = ( — 

x T7.7 



vl/3 



is the Sedov length, A 



(16) 
cST is the 



width of the shell (ST is the duration of GRB) and To is 
the initial Lorentz factor of the ejecta. If £ < 1, the reverse 
shock becomes relativistic before it crosses the shell; oth- 
erwise (£ > 1), the reverse shock remains Newtonian or at 
best mildly relativistic during the whole energy extraction 
process. For GRB990123, we find £ ~ 0.3 < 1. So the re- 



verse shock of GRB990123 had become relativistic before it 
crossed the shell, consistent with our calculated result. 



4 CONCLUSIONS AND DISCUSSIONS 

We have constrained some intrinsic parameters, such as the 
magnetic energy density fraction (es), the electron energy 
density fraction (e e ), the isotropic energy in the adiabatic 
forward shock £52 and the ambient density n. Due to the 
lack of the value of the self-absorption frequency v a , we made 
an assumption that e e of GRB990123 is the same as that of 
GRB970508, then astonishingly find that the inferred value 
of es is also nearly equal to that of GRB970508. This result 
favours the argument proposed by WG99 that the magnetic 
energy fraction and the electron density fraction may be 
universal parameters. 

Another two important intrinsic parameters of 
GRB990123 are also inferred from the optical flash informa- 
tion: the initial Lorentz factor Fo and the Lorentz factor Y a 
at the prompt optical emission peak time of the ejecta. They 
are: To = 1200, Y a = 300. Our inferred value of the To is 
six times larger than that obtained by Sari & Piran (1999b), 
who used the ambient density n inferred for GRB970508. A 
larger initial Lorentz factor is reasonable in consideration of 
the huge energy of this burst. The Lorentz factor of the re- 
verse shock at the optical flash peak time is Y rs ~ £&■ ~ 4, 
which shows that the reverse shock had become relativistic 
rather than mildly relativistic before it crossed the entire 
ejecta shell. This result is in agreement with the criterion 
presented by Sari & Piran (1995) to judge the RRS case or 
NRS case. 

Prompt optical flash has added another dimension to 
GRB astronomy. Prompt observations in the optical band 
during and immediately after GRB may provide more and 
more events of optical flash in the near future, and they 
will enable us to make more detailed analyses, make more 
precise determination of intrinsic parameters and test the 
reverse — forward external shock model. 
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